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Photoinduced dynamics of charge density and lattice displacements is calculated by solving 
the time-dependent Schrodinger equation for a one-dimensional extended Peierls-Hubbard 
model with alternating potentials for the mixed-stack organic charge-transfer complex, TTF- 
CA. A pulse of oscillating electric field is incorporated into the Peierls phase of the transfer 
integral. The frequency, the amplitude, and the duration of the pulse are varied to study the 
nonlinear and cooperative character of the photoinduced transition. When the dimcrizcd ionic 
phase is photoexcited, the threshold behavior is clearly observed by plotting the final ionicity 
as a function of the increment of the total energy. Above the threshold photoexcitation, the 
electronic state reaches the neutral one with equidistant molecules after the electric field is 
turned off. The transition is initiated by nucleation of a metastable neutral domain, for which 
an electric field with frequency below the linear absorption peak is more effective than that 
at the peak. When the pulse is strong and short, the charge transfer takes place on the same 
time scale with the disappearance of dimerization. As the pulse becomes weak and long, the 
dimerization-induced polarization is disordered to restore the inversion symmetry on average 
before the charge transfer takes place to bring the system neutral. Thus, a paraelectric ionic 
phase is transiently realized by a weak electric field. It is shown that infrared light also 
induces the ionic-to-neutral transition, which is characterized by the threshold behavior. 

KEYWORDS: photoinduced phase transition, neutral-ionic transition, TTF-CA, charge-transfer 
complex, time-dependent Schrodinger equation, dynamics 

1. Introduction 

Photoinduced phase transitions^^ are intriguing phenomena, where relatively weak pertur- 
bation causes a macroscopic change in magnetic, dielectric, optical and structural properties. 
It would be important to study their mechanisms to explore possibilities of dynamically con- 
trolling these properties by making use of the cooperativity that often allows spontaneous 
symmetry breaking. There are also possibilities of separately controlling some of these elec- 
tronic properties on the respective time scales. 

Now many materials are known to show photoinduced phase transitions. The characters 
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of the transitions arc very different, depending on the relative importance of the interaction 
responsible for the cooperativity to couplings with environmental surroundings. In the spin 
crossover complexes, the spin state changes very slowly compared with molecular vibrations 
coupled with it. The initial and final spin states have their respective diabatic potentials as 
a function of the so-called interaction mode, whose crossing gives the barrier height to be 
overcome by thermal lattice fluctuations. The details of lattice dynamics is unimportant, 
so that we can safely integrate over the thermally fluctuating lattice displacements.^) The 
evolution of the spin state is thus described in a stochastic manner, e.g., by using a master 
equation^' or Monte Carlo simulations.®' ^-^ In these materials, the direct interaction among 
rather distant d-electron spins is much weaker than the effective interaction, which is mediated 
by the surrounding lattice displacements. 

In contrast, in the mixed-stack organic charge-transfer complex, TTF-CA 
(TTF=tetrathiafulvalene, CA=chloranil),^'^"-'^'^) whose model is studied in this paper, 
the dynamics of charge density and that of dimerization- induced, staggered lattice displace- 
ments are strongly coupled. They are governed by the electron-lattice and electron-electron 
interactions^^-* that are responsible for the electronic and structural properties in equilibrium. 
The electronic state substantially evolves according to these interactions before the energy 
supplied by the photoexcitation significantly dissipates into unimportant degrees of freedom 
such as intramolecular vibrations. In other words, on the time scale of observation, which is 
rather short for TTF-CA, memories are not easily lost as in the slowly evolving spin crossover 
complexes. The rapidncss of the electronic evolution is generally expected for a low energy 
barrier, i.e., when the hysteresis is small or when the state is near the instability within the 
hysteresis. 

Indeed the coherent motion of the macroscopic neutral-ionic domain boundary is ob- 
served, which would disappear if the energy dissipation were significant. The evolution of 
the electronic state should thus be described in a deterministic manner, e.g., by combining 
the time-dependent Schrodinger equation for the electronic wave function and the classical 
equation of motion for the lattice displacements. We have actually described in ref. 18 the 
charge-lattice dynamics that does not follow the adiabatic potential that is projected on the 
plane spanned by the shape and size of a metastable domain.^^^ In this paper, we again use 
the Hartree-Fock (HF) approximation and introduce photoexcitations not by simply changing 
the occupancy of the HF orbitals^^^ but by incorporating a pulse of oscillating electric field 
into the Peierls phase of the transfer integral. Thus, the frequency, the amplitude, and the 
duration of the pulse can be independently varied. 

Concerning the issue of whether the evolution is stochastic or deterministic, we note that 
the photoinduced charge-density-wave(CDW)-to-charge-polarization(CP) phase transition in 
one-dimensional halogen-bridged binuclear platinum complexes show both the stochastic and 
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deterministic characters. ^'^^ The dependence of the threshold photoexcitation on the electron- 
lattice coupling is intuitively explained by the dependence of the barrier height. ^"^^ Transition 
processes between the two electronic states may be described in a stochastic manner with the 
help of their diabatic potentials if only this fact is to be explained. However, the difficulty of the 
CP-to-CDW phase transition compared with the CDW-to-CP one needs detail explanations 
based on charge-transfer processes in the respective electronic states. These processes are 
easily missed if one uses a stochastic approach. 

Among the materials that show rather rapid photoinduced charge-lattice-coupled dynam- 
ics, the TTF-CA complex would be the most extensively studied, ^'^"^^^ which allows compar- 
isons with model calculations. In this paper, we discuss the transition from the ionic phase to 
the neutral one by intra-chain charge-transfer photoexcitations, for which the coherent mo- 
tion of the macroscopic neutral-ionic domain boundary is observed as well as the threshold 
photoexcitation density. -"^^^ In the case of intramolecular excitations, a new ionic phase with 
disordered polarizations is suggested to appear by comparing the photoreflectance and the 
second-harmonic-generation (SHG) signal^^^ and by the x-ray diffraction.^''^ The SHG inten- 
sity, which corresponds to the degree of the inversion-symmetry breaking, decreases after the 
photoirradiation much faster than the ionicity estimated from the photoreflectance. The x-ray 
diffraction gives more direct evidence, but the error bar is not so small at present. 

In this paper, we extend the previous study, where the dynamics of charge density 
and that of lattice displacements during the photoinduced ionic-to-neutral phase transition is 
discussed in a model for TTF-CA. Here we directly incorporate a pulse of oscillating electric 
field. Thus, we can now study the dependence on the frequency, the amplitude, and the 
duration of the pulse. Furthermore, this approach allows us to discuss a possibility for an 
infrared-light- induced transition without producing excitons, as shown in the present paper. 
It also allows us to discuss the effect of a double pulse with different intervals, as demonstrated 
later in paper III. 

2. Extended Peierls-Hubbard Model with Alternating Potentials 

The highest occupied molecular orbital (HOMO) at the donor site and the lowest unoc- 
cupied molecular orbital (LUMO) at the acceptor site are necessary to describe the electronic 
state. In the limit of vanishing transfer integral to ~^ 0, the HOMO is doubly occupied in the 
neutral phase, while both of the HOMO and the LUMO are singly occupied in the ionic phase. 
Different types of electron-lattice couplings have so far been employed to explain the dimer- 
ization.^^"^'*) Here we assume that the Coulomb interaction strength is modified by the lattice 
displacement.^^) We use a one-dimensional extended Peierls-Hubbard model with alternating 
potentials at half filling, which is a single-chain analogue of the model used in ref. 19, 

H = Hel + i/lat , (1) 
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where, cj^ ici,a) is the creation (annihilation) operator of a 7r-electron with spin a at site I, 
f^i,(T = cJct'^^>'^' ~ "',T "I" i^ dimensionless lattice displacement of the Zth molecule 

along the chain from its equidistant position, and yi = u^+i — The distance between the 
/th and {I + l)th molecules is then given by ri = ro(l + it/+i — it/), where tq is the averaged 
distance between the neighboring molecules along the chain. The donor and acceptor molecules 
are located at the odd and even sites, respectively. The nearest-neighbor repulsion strength 
between the Zth and {I + l)th sites Vi depends on the bond length yi, Vi = V + f32yfi where V 
is for the regular lattice, and P2 is the quadratic coefficient. The linear coefficient is neglected 
for simplicity. The parameter Iq denotes the nearest-neighbor transfer integral, U the on-site 
repulsion strength, and d the level difference between the HOMO and the LUMO in the neutral 
phase with to = 0- The elastic energy is expanded up to the fourth order: the parameters ki 
and ^2 are the linear and nonlinear elastic constants. The mass of the Zth molecule is denoted 
by nil. The number of sites is denoted by N. 

In the unrestricted Hartrec-Fock (HF) approximation, we iteratively solve the eigenvalue 
equation to obtain the ground state that is self-consistent with the one-body densities {c^ ^ci^a) 
and {cl ^ci^i^cr) and the lattice displacements yi, which satisfy the Hellmann-Feynman theorem. 
We then add random numbers to the initial yi and iii values according to the Boltzmann 
distribution at a fictitious temperature T. With thus modified lattice displacements, the HF 
Hamiltonian is diagonalized once again by imposing the self-consistency on the one-body 
densities only. 

Photoexcitations are introduced by modifying the transfer integral in the kinetic part of 
the electronic Hamiltonian, 

N 

- YI (io(t)4<,Q+i,, + h.c.) , (4) 



with the Peierls phase. 
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where e is the absolute value of the electronic charge, and c the light velocity. The time- 
dependent vector potential A{t) is related to the electric field E{t) by 



with amplitude ^ext and frequency cjext for < t < -/Vext^ext with integer A?ext- E{t) is zero 
otherwise. The period is denoted by Text, which is given by Text = 27r/a;ext- 

Now we mention how the collective nature is incorporated into this photoexcitation. The 
excitonic effect, namely, the attraction between an excited electron and a hole to form a bound 
pair can be included, e.g., in the random phase approximation (RPA). The RPA treats all 
the linear quantum fluctuations around the unrestricted HF solution.^^^ These fluctuations 
correspond to the infinitesimal deviations from the unrestricted HF state, all of which as well 
as their interactions are included in the time-dependent HF theory. Deviations from the un- 
restricted HF state are no longer infinitesimal when we solve the time-dependent Schrodinger 
equation after the photoexcitation. Thus, the excitations generally have a finite lifetime in 
the present numerical calculations. The present treatment of photoexcitation does include the 
excitonic effect from the beginning, as clearly shown later in the linear absorption spectrum. 

For the electronic part, we thus solve the time-dependent Schrodinger equation, as ex- 
plained in ref. 18, during and after the photoexcitation. The exponential operator is decom- 
posed,^^) in such a way that the decomposition is accurate to the order of {At)^?"^^ For the 
lattice part, we solve the classical equation of motion, as explained in ref. 18. The leapfrog 
method is employed, which is accurate to the order of (At)^. After the electric field is turned 
off, the total energy is conserved within the numerical accuracy. 

3. Results and Discussions 

We use A^=100, to=0.17eV, ?7=1.528eV, y=0.604eV, d=2.716eV, /32=8.54eV, A;i=4.86eV, 
A;2=3400eV, and the bare phonon energy Wopt = (l/ro)(2A;i/mr.)-'^/^=0.0192eV, and impose the 
periodic boundary condition. The reduced mass rrir is defined as rrir = mBmA/{mY) + jtza) 
with mD for the donor molecule and tua for the acceptor molecule. With these parameters, 
the dimerized ionic phase is stable and the neutral phase is metastable. These parameters 
are the same as in ref. 18, so that the bare phonon energy c^opt used here is about five times 
higher than the optical phonon energy of the TTF-CA complex. The results are very similar 
for smaller Uopt values including the one for the TTF-CA complex,^^) so that the conclusion 
is not altered by specific values of Wopt- Simply because it reduces the computational burden, 
we adopt this value. The ionicity is defined as p = 1 -|- (1/-/V) I)'('T'z)- The staggered 

lattice displacement is defined as yst = (V-^) 




(6) 



where we consider a pulse field. 
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Fig. 1. Linear absorption, as a function of Woxt/'^opt in the excitonic range, in the ionic phase with 
lattice fluctuations at T/t=lQ~^. The applied electric field is very weak, ea-Eext/wopt=0.01, and 
lasts long, Afext=1000. 



3.1 Linear absorption 

First, we turn on a very weak electric field for a long time and observe the increment of the 
total energy. It is plotted in Fig. 1, as a function of the frequency of the electric field. Since we 
here use a;opt=0.0192eV, the position of the absorption peak is Wext — 32a;opt — 0.61eV, which 
corresponds to the peak in the imaginary part of the current-current correlation function 
obtained by the RPA.^^) The calculation of the linear response in this subsection is easily 
extended to those of nonlinear responses including the electromodulation spectrum.^^) Initial 
lattice fluctuations at r/f=10~^, i.e., T/uopt ::i^0.0089 broaden the absorption spectrum. The 
spectrum is smooth especially when we apply a very weak field for a very long time, which 
would effectively average random distributions of the lattice displacements. It is reminiscent of 
the ergodic theorem that guarantees the replacement of the distribution average by the time 
average. Since the ionic phase is directly affected by the nearest-neighbor repulsion Vi and 
thus by the lattice displacements, its absorption spectrum is much more sensitive to lattice 
fluctuations than that in the neutral phase in the present model. 

3.2 Threshold behavior 

To see how the linear and nonlinear absorption spectra are different, we now use a stronger 
and shorter pulse. By adopting ea£'ext/'*^opt=5 and iVext=20, we again calculate the increment 
of the total energy. In Fig. 2, the increment is divided by the frequency of the electric field, 
which corresponds to the number of absorbed photons in the 100-site chain. The strong 
absorption is extended to rather low frequencies [Fig. 2(a)]. In the frequency range where 
many photons are absorbed, the system drastically changes, namely, the transition takes 
place to the neutral phase. Accordingly, the ionicity finally drops [Fig. 2(b)]. 
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Fig. 2. (a) Number of absorbed photons and (b) final ionicity, as a function of Wext/^^opt in the 
excitonic range, when the electric field with ea£'ext/'^opt=5 and -A/ext=20 is applied to the ionic 
phase at T/t=10-^. 



By setting the frequency of the electric field a;ext near the linear-absorption peak, 
ujext/^opt = 32, we vary the field strength (eaE'ext)^ and the pulse duration NextText in Fig- 3 
to see how they change the number of absorbed photons and the final ionicity. When the 
pulse is weak or short, the number of absorbed photons is almost proportional to the field 
strength or the pulse duration. This situation is smoothly connected to the linear absorption 
shown in Fig. 1 (i.e., where the field is infinitesimally small). When the number of photons 
is more than about eleven, the absorption strongly deviates from the linearity. With stronger 
fields, the absorption tends to be saturated [Fig. 3(a)]. Then, the final ionicity quickly drops, 
demonstrating that the transition to the neutral phase takes place above the threshold ab- 
sorption of about eleven photons in the 100-site chain [Fig. 3(b)]. In the previous study,^^^ 
photoexcitations were introduced by simply exchanging the occupancy of the initial HF or- 
bitals around the Fermi level. With the same random numbers added to the lattice variables 
(i.e., in the same initial condition), the threshold absorption was about four photons in the 
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Fig. 3. (a) Number of absorbed photons as a function of the field strength (eaE'ext/'^opt)^, and (b) 
the corresponding final ionicity, for A^ext=2- (c) Number of absorbed photons as a function of 
the pulse duration Next, and (d) the corresponding final ionicity, for ea£'ext/wopt=7. The electric 
field of frequency a;ext/wopt=32 near the linear-absorption peak is applied to the ionic phase at 
T/t=10-^. 



100-site chain. The more realistic, oscillating electric field is therefore less effective than the 
simple exchange of the orbital occupancy. The absorption increases with the pulse duration 
more rapidly when the number of absorbed photons is more than about eleven [Fig. 3(c)]. 
The threshold absorption is again given by about eleven photons in the 100-site chain. Above 
the threshold pulse duration, the final ionicity gradually increases [Fig. 3(d)]. This indicates 
that, after the transition to the neutral phase is completed, further application of the oscil- 
lating electric field leads to further absorption of photons [Fig. 3(c)], which results in random 
distribution of electrons in the present calculations without dissipation. Because the orbitals 
are delocalized, the completely random distribution of electrons in the thermodynamic limit 
gives one electron per site, corresponding to p=l, at half filling. In all cases, such a disordered 
state with p ~1 is always reached by sufficiently long application of the oscillating electric 
field. 
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4. (a) Number of absorbed photons as a function of the field strength (eaE'ext/'^opt)^, and (b) 
the corresponding final ionicity, for A^ext=3- (c) Number of absorbed photons as a function of 
the pulse duration A'ext, and (d) the corresponding final ionicity, for ea£'ext/'^opt=5. The electric 
field of frequency a;ext/'^opt=28 below the linear- absorption peak is applied to the ionic phase at 
T/t=lQ-^. 
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Now we use a frequency Wext below the linear-absorption peak, Wext/^'^opt = 28, in Fig. 4 to 
see how the above behavior is modified when the field is off resonance. It is generally similar 
to that when the field is on resonance. Namely, when the number of photons is more than 
about eleven, the absorption strongly deviates from the linearity [Figs. 4(a) and 4(c)], and 
the transition to the neutral phase takes place above this threshold absorption [Figs. 4(b) 
and 4(d)]. However, to go beyond the threshold absorption, the pulse duration with A?ext=4 
is sufficient for ea^ext/^^opt=5 [Figs. 4(c) and 4(d)], which is compared with A?ext=25 for 
ea£'ext/'^opt=7 [Figs. 3(c) and 3(d)] on resonance. The frequency a;ext/'^opt=28 is located 
in the low-frequency tail of the linear absorption. The efficiency in the off-resonant condition 
comes from the fact that excitons localized due to lattice fluctuations nucleate neutral domains 
more easily than delocalized excitons. This fact will be discussed later again. 

At the threshold absorption, we estimate the ratio of the energy needed for the phase 
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transition to the energy supplied by the photoirradiation. The energy difference between the 
stable ionic phase and the metastable neutral phase is 0.16eV for the 100-site chain (Fig. 1 
of ref. 18). The energy to create two neutral- ionic domain walls, which are accompanied 
with each metastable domain, is about 0.3eV according to Fig. 10 of ref. 19. (It is about 
a half of that in the strong-coupling limit, y=0.604eV.^'^)) Meanwhile, the frequency of the 
oscillating electric field used here is either 32LJopt — 0.61eV, which corresponds to the energy 
of a delocalized exciton at the linear absorption peak, or 28a;opt — 0.54eV. Therefore, a single 
photon is energetically sufficient to bring about the transition. The increment of the total 
energy due to initially given lattice fluctuations is 0.07eV before the photoirradiation, so that 
it is much smaller than the values above. Just above the threshold absorption, the lattice 
kinetic energy is finally in the range between O.leV and 0.2eV for the 100-site chains, and the 
lattice elastic energy is comparable with it, in all the cases shown in Figs. 3, 4 and 11. 

These values indicate that more than 90% of the energy supplied by the photoirradiation is 
absorbed into the electronic state, leading to deviation from the static self-consistent solution 
in the metastable neutral phase. The final ionicity apparently deviates from 0.2 in the static 
solution, as shown in Fig. 10 of ref. 18. Thus, the threshold density of absorbed photons (eleven 
photons in the present case, four photons by the simple exchange of the orbital occupancy,^®) 
in the 100-site chain) is not important. Important is the fact that less than 10% (about 
20% if the orbital occupancy is exchanged) of the supplied energy is used to bring about the 
transition along the fact that this ratio is insensitive to whether the field is on resonance or off 
resonance. Later, we will show that the ratio of the energy difference between the two static 
electronic states to the energy actually supplied by the photoirradiation somewhat depends 
on the field strength or the pulse duration. The ratio becomes low as the field is weakened. 
In general, as the energy dissipates more into processes other than charge transfers, the ratio 
becomes low. To theoretically estimate the ratio in this many-body system is far beyond the 
scope of the present study. It will need accumulations of numerical calculations. 

We have calculated the evolution of the ionicity with various strengths and durations of 
pulses to find their relation at the threshold. The pulse duration A'ext just above and below 
the threshold for the transition is plotted in Fig. 5(a), as a function of the field strength 
(ea£^ext/'i^opt)^- In general, to describe the dynamics of the density of excited states under 
photoirradiation, a master equation is useful when thermal transitions are so dominant that 
the evolution is described in a stochastic manner.^) In such a case, it is known that the time 
required for a transition icross is proportional to (/ — Ith)~^^'^ when the field strength / is just 
above the threshold /th-^^ Then, we also plot N~J^ in Fig. 5(b), as a function of the strength. If 
the same law held as in the stochastic case, would be a linear function of (ea£^ext/'^opt)^ 
and would vanish at a finite value of (ea-Ecxt/i^^opt)^; for which the transition would never take 
place no matter how long the field is applied. However, the present result shows a nonlinear 
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curve, suggesting the absence of ith although it is very hard to demonstrate it numerically. 
Even if the field is very weak, e.g., (ea£^ext/'*^opt)^=l) the transition finally takes place after 
the field is applied for a very long time. Thus, in the deterministic case, the conversion occurs 
more efficiently than in the stochastic case. In the experimental condition, if a weak field is 
applied, a portion of energy would flow into the environment (i.e., outside the model system) 
before the transition occurs, so that the efficiency would become lower. As the exchange of 
energy between the system and the environment dominates the direct exchange inside the 
system, the efficiency would approach the value predicted in the stochastic case. 

The evolution of the ionicity is shown in Fig. 6(a) with different field strengths 
{eaEext /(^opt)^ , and in Fig. 6(b) with different pulse durations Next- Note that the oscillating 
electric field is applied in Fig. 6(a) for < Wopt* < -/Vextt^optrext=-/Vext x 27r/[a;ext/t^opt]=3 x 
27r/28=0.67. With weak fields below the threshold, the absorbed energy partially goes into the 
lattice oscillations, sometimes reversing the polarization of some ionic domains, and partially 
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into the electronic part, increasing the deviation from the self-consistency condition, but the 
charge is not transferred so much between the donor molecules and the acceptor molecules. 
Above the threshold strength, a neutral domain is locally created before the field is turned 
off at t=NcxtTcxt- At this time, the spatially averaged ionicity is only slightly reduced if the 
field is just above the threshold. After that, the neutral domain grows, although the growth 
is not always monotonous due to the competition with surrounding ionic domains. In rare 
cases, as discussed later, an initially created neutral domain is completely annihilated and 
another neutral domain is created in a different place, which grows and finally covers the 
whole system. With stronger fields much above the threshold, the initially created neutral 
domain monotonously grows and considerably changes the spatially averaged ionicity already 
at t=NcxtTcxt- Similar evolution is observed by increasing the pulse duration. With short du- 
rations below the threshold, the ionicity is almost unchanged. Above the threshold duration, 
a neutral domain is created and grows after the field is turned off. In Fig. 6(b), the curves for 
-^cxt=4 and 6 show this situation. In both cases, the neutral domain is transiently suppressed 
by surrounding ionic domains after t=NcxtTcxt- With a longer duration (A'cxt=6), the system 
absorbs more energy and the transient suppression of the neutral domain is weaker. In any 
case, the threshold behavior shown in Figs. 3 and 4 is a consequence of the fact that the 
ionicity either reaches a value in the neutral phase or almost returns to the initial value (i.e., 
no intermediate values between them) after the field is turned off, as clearly shown in Fig. 6. 

The space and time evolution of the ionicity and the staggered lattice displacement is 
shown in Fig. 7. The horizontal component of each bar represents the local ionicity pi defined 
as pi = 1 + {—ly {— (ni^i) + 2{ni) — (n;+i))/4. The vertical component gives the local staggered 
lattice displacement tjst i defined as i/st i = (— 1)' (—?/;_ i -|- 2yi — ?/;+i)/4. The bars arc shown on 
all sites and selected times. Now the oscillating electric field is applied for < LOoptt < 33.7. 
The initial ionic phase at t=0 has uniformly downward bars (not shown). Note that, for 
"^opt^ < 33.7, some ionic domains have reversed polarizations. When the field is turned off, the 
local ionicity is weakened around 90 < / < 100. This local area can be regarded as the nucleus 
or the seed of a neutral domain, although the ionicity of this area at o^opt^ = 33.7 is still larger 
than in the neutral domain appearing later. Once the neutral domain appears here, it quickly 
grows and covers the 100-site system within one or two periods of the optical lattice vibration 
(whose period is about 4/u;opt)- The boundary between the ionic and neutral domains is 
clearly seen because the ionicity changes discontinuously at the boundary. Immediately after 
the transition is completed, the ionicity is low but still spatially fluctuates. With increasing 
Wopt^) the ionicity becomes spatially uniform. 

In order to analyze the time-dependent electronic state, we obtained Fourier transforms 
of the ionicity in ref. 18 to compare with them the current excitation spectra of the static 
neutral state in the RPA. The RPA is equivalent to the time-dependent HF approximation 
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Fig. 6. Time dependence of the ionicity, (a) with different field strengths (eaE'ext/'^opt)^ for A^ext=3, 
and (b) with different pulse durations A/ext for eai^ext/'^opt=5- The electric field of frequency 
t^ext/t^opt=28 below the linear-absorption peak is applied [(a) for < Wopt* < -^ext'^opt^ext=0.67, 
and (b) for < Wopt* < -/Vexti^optT'ext=0.45, 0.67, 0.90 and 1.35 with A^oxt=2, 3, 4 and 6, respec- 
tively] to the ionic phase at T/t=10~^. The final state is either ionic or neutral. 



used here but limited to infinitesimal deviations from the static HF solution, so that the linear 
absorption peak in Fig. 1 corresponds to the peak in the RPA spectrum due to delocalized 
excitons. After a sufficiently large number of electrons are excited, the time-dependent ionicity 
has high-frequency components due to delocalized excitons and low-frequency components due 
to complex motion of the neutral-ionic domain walls. When the number of excited electrons 
decreases, the correspondence to delocalized excitons becomes worse because excitons are more 
sensitively affected by lattice fluctuations.^^) Since the frequency of the electric field applied in 
Fig. 7 is below the peak in Fig. 1, the photoirradiation produces excitons localized by lattice 
fluctuations. 

Although the time-dependent Schrodinger equation governs the deterministic time- 
evolution, the positions where excitons are produced are distributed in a stochastic manner 
according to the random numbers initially added to the lattice variables. As seeds of stable 
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Fig. 7. Correlation between the staggered lattice displacement yst i (the vertical component of the 
bar) and the ionicity pi (the horizontal component of the bar), as a function of the site index 
I and the elapsing time t multiplied by uoopt- The electric field with ea_Bcxt/t^opt=2, -/Vcxt = 150 
and of frequency a;cxt/wopt=28 below the linear- absorption peak is applied (for < Wopt^ < 
-^ext't^optText=33.7) to the ionic phase at T/t=lQ~^. The final state is neutral. 



domains generally trigger a first-order transition, localized excitons promote the transition. 
Indeed, the localized excitons more efficiently nucleate neutral domains and trigger the tran- 
sition (Fig. 4) than the delocalized excitons (Fig. 3). Neutral domains are created by the 
localized excitons, so that they are also distributed at random. Although these nucleation 
processes are stochastic, the neutral domains grow in a deterministic manner as long as the 
neutral-ionic domain walls do not collide with ionic-anti-ionic (I-I) solitons, which are bound- 
aries between different polarizations inside the ionic domain, as discussed in ref. 18. That 
is why the coherence can appear, as discussed later in paper III. Namely, if the stochastic 
processes dominated the time-evolution of the system, the coherence would be lost. 

3. 3 Strong vs. weak pulses 

To see how the electronic dynamics and the lattice dynamics are correlated and how 
their correlation can be controlled by changing the parameters of the oscillating electric field, 
we use different strengths and durations of pulses to calculate the time evolution of the 
ionicity and of the staggered lattice displacement. With a strong and short pulse, the electronic 
dynamics and the lattice dynamics are closely correlated. One example is shown in Fig. 8. On a 
magnified scale, one sees small-amplitude and very rapid modulation of the ionicity [Fig. 8(a)] 
related with excitonic excitations, in contrast to the slow oscillation of the staggered lattice 
displacement [Fig. 8(b)]. However, the overall evolutions of the ionicity and of the staggered 
lattice displacement are similar to each other. 

As the pulse becomes weaker and longer, their evolutions become different from each other, 
as shown in Fig. 9. The ionicity initially drops a little bit, but it remains large for a while 
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Fig. 8. Time dependence of (a) the ionicity and (b) the staggered lattice displacement when the pulse 
of frequency a;ext/'^opt=28 is rather strong, eai^ext/'^opt=3, and short, A^ext=10 (i.e., < Wopt* < 
-^extWoptT'ext=2.24). The initial state is ionic at T/t=10~^. The ionicity and the staggered lattice 
displacement change on the same time scale to reach the neutral phase. 

[Fig. 9(a)]. During the same period, the staggered lattice displacement steadily decreases and 
almost vanishes [Fig. 9(b)]: the polarization of some ionic domain is reversed, then that of 
another domain is reversed, and so on. In other words, I-I solitons separating ionic domains 
with different polarizations are created in a random manner during this period. When the 
staggered lattice displacement vanishes on average, the system is regarded as in the para- 
electric ionic phase. It should be noted that this phase is known to appear in TTF-CA only 
under high pressure^^^ in thermal equilibrium. The mechanism of its appearance is theoreti- 
cally clarified in the context of the electrostriction effect.^^^ Furthermore, after intra-molecular 
photoexcitations, this phase is suggested to appear by comparing the photorefiectance and the 
second-harmonic-generation signal^^^ and by the x-ray diffraction.^^) The present numerical 
result suggests that, even by intra-chain charge-transfer excitations, it is possible to tran- 
siently create the paraelectric ionic phase if the pulse is so weak and long that I-I solitons are 
created before the charge transfer really occurs. When the pulse is weak, the supplied energy 
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Fig. 9. Time dependence of (a) the ionicity and (b) the staggered lattice displacement when the pulse 
of frequency a;ext/'^opt=28 is weaker, ea£'ext/'^opt=2, and longer, Ne^=l^Q (i.e., < Wopt* < 
-^ext'^opt^ext=33.7), than in Fig. 8. The initial state is ionic at T/t=lQ~^. The staggered lattice 
displacement decays faster than the ionicity to reach the neutral phase. 

is not used to directly transfer charge density along the chain. By intramolecular excitations, 
the supplied energy would not be directly transferred to intra-chain charge transfer processes, 
either. In this sense, the paraelectric ionic phase may appear by a common mechanism. If the 
electric field continues to be applied, this phase is finally converted into the neutral phase. 

The weak-pulse case where the supplied energy dissipates more into I-I solitons is similar 
to the large-T case in ref. 18. When the initial lattice fluctuations are small, few I-I solitons 
are present. The neutral domain grows smoothly and rapidly (Fig. 5 of ref. 18). Otherwise 
many I-I solitons collide with the neutral-ionic domain walls to make the motion of the latter 
irregular. Then, the growth of the neutral domain becomes rather chaotic (Fig. 7 of ref. 18). 
In this paper, many I-I solitons are created when the weak field is applied for a long time, 
and the phases of the staggered lattice displacements (i.e., the polarizations) are disordered 
to average out the broken inversion symmetry. Accordingly, the growth of the neutral domain 
becomes rather irregular, e.g., in the case of Fig. 10, although it is invisible on this scale. 
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Fig. 10. Time dependence of (a) the ionicity and (b) the staggered lattice displacement when the 
pulse of frequency a;ext/'^opt=28 is much weaker, ea£'ext/'^opt=l) and much longer, Afext=6800 
(i.e., < cOoptt < A?extt^opt^ext=1530), than in Fig. 8. The initial state is ionic at T/t=10~^. 
The ionic phase without long-ranged dimerization pattern (thus paraelectric) lasts long before the 
neutral phase is reached. 



When the very weak electric field is applied continuously, the ferroelectric ionic phase 
survives for a while, as shown in Fig. 10. Then, the ionicity drops a little bit, but further 
charge transfer does not take place for a long time [Fig. 10(a)], while the absorbed energy 
slowly disorders the polarization of the ionic phase to reduce and finally destroy the long- 
range order of the staggered lattice displacement [Fig. 10(b)]. It takes a long time to suppress 
the staggered lattice displacement and an even longer time to create a neutral domain. The 
paraelectric ionic phase has a long lifetime. It can be made infinite within the computational 
time scales if the field is turned off before the neutral domain is created. Once the neutral 
domain is created, it grows as a whole. The I-I solitons obstruct the growth of the neutral 
domain,^^) but they cannot completely stop it. 

As stated above, when the weak field is applied for long, the supplied energy dissipates 
more into I-I solitons and others, so that it is not used to directly transfer charge density 
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Fig. 11. Final ionicity as a function of the number of absorbed photons, (a) for pulses of intermediate 
strength, eai^ext/'^opt=3, and (b) for weak pulses, ea£'ext/'^opt=l and 2. The electric field of 
frequency a;ext/t^opt=28 is applied to the ionic phase at T/t=10~^. 



along the chain. It somewhat modifies the threshold relation. The final ionicity is then plotted 
in Fig. 11, as a function of the number of absorbed photons, for intermediate to weak pulses. 
The threshold absorption is indeed shifted to larger values. Namely, it is about 11 photons 
per 100 sites for ea£'ext/<^opt=5 [Fig. 4(d)], 11 or 13 photons for ea£'ext/'^opt=3 [Fig. 11(a)], 
15 photons for ea£^ext/'*^opt=2, and 27 photons for ea£'ext/'^opt=l [Fig- 11(b)]. In Fig. 11(a), 
the final ionicity is not a monotonous function. As the pulse duration increases, it once drops 
and returns to a large value before it finally drops again to a small value. This is because 
the initially created neutral domain is annihilated by continuing photoirradiation. It never 
happens if the field is strong. For intermediate to weak pulses, larger dissipation or larger 
excess energy can disturb the growth of the neutral domain. The second neutral domain, 
which finally grows and covers the whole system, is created in a different place from the first 
one. 

Here we estimate the energy flow again. It is true that many I-I solitons are produced 
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for weak and long pulses, but their creation energies are much smaller than those of the 
neutral-ionic domain walls mentioned before. Both of the lattice kinetic and elastic energies 
are finally in the range between O.leV and 0.2eV for the 100-site chain. It implies that most 
of the supplied energy is absorbed into the electronic state. Consequently, the final ionicity 
deviates from 0.20 for the static self-consistent solution in the metastable neutral phase: it is 
0.32 for ea£'ext/'^opt=3 (Fig. 8), 0.35 for ea£'ext/'^opt=2 (Fig. 9), and 0.43 for ea£'ext/wopt=l 
(Fig. 10). 

3.4 Infrared pulses 

Because of the electron-lattice coupling, photons whose frequencies arc about that of the 
optical lattice vibration are also absorbed and can contribute to the phase transition. The 
transition to the neutral phase is induced by sufficiently strong irradiation of the infrared 
light. By setting the frequency of the electric field ojcxt comparable with that of the optical 
lattice vibration, we vary the field strength (eaEcxt)'^ and the pulse duration A'^cxt^'cxt- Below 
threshold strength or duration, photons arc hardly absorbed. Above the threshold values, 
many photons are absorbed. The important point here is again the fact that, once substantial 
absorption takes place, the ionicity substantially decreases. Besides the fact that the relation 
between the number of absorbed photons and either the field strength or the pulse duration 
is now strongly nonlinear (not shown), the presence of the threshold and the behavior above 
the threshold are very similar to those in the case of excitonic excitations. 

4. Conclusions 

In the quasi-one-dimensional mixed-stack organic charge-transfer complex, TTF-CA, pho- 
toirradiation is known to trigger transitions from the ionic to neutral phases and from the 
neutral to ionic phases. In the former transition, the coherent motion of the macroscopic 
neutral-ionic domain boundary is observed as well as the threshold density of absorbed 
photons. Coherence will be discussed later. To deterministically approach the dynamics 
of charge density coupled with lattice displacements during this transition, we employ the 
one-dimensional extended Peierls-Hubbard model with alternating potentials and numerically 
solve the time-dependent Schrodinger equation for the electronic wave function at the mean- 
fled level combined with the classical equation of motion for the lattice displacements. To 
treat the photoexcitations more realistically than in the previous study, we incorporate a 
pulse of oscillating electric field into the Peierls phase of the transfer integral and vary the 
frequency, the amplitude, and the duration of the pulse independently. 

When the dimerized ionic phase is photoexcited, we reproduce the threshold behavior 
in the final ionicity as a function of the increment of the total energy, which is regarded as 
proportional to the density of absorbed photons. The threshold photon density rather weakly 
depends on the amplitude and the duration of the pulse. However, when the final ionicity 
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is plotted as a function of the amplitude or duration of the pulse, we found that the pulse 
with frequency below the linear absorption peak is more effective than that at the peak. The 
phase transition proceeds by nucleation of a metastablc neutral domain and its rapid growth. 
As a consequence, after the electric field is turned off, the electronic state either reaches the 
neutral one (above the threshold photocxcitation) or returns to the ionic one with excess 
energy supplied by photoexcitations (below the threshold), but basically no intermediate 
state is produced. The efficiency of the pulse with frequency below the linear absorption peak 
originates from the fact that lattice-fluctuation-induced localized excitons are more effective 
on nucleating the neutral domain than delocalized excitons. 

Here we discuss coherence. Below and in the following papers, we use this word in the 
numerical context to describe a state where the phases of the staggered lattice displacements 
are almost uniform in a wide spatial and temporal region inside the ionic domain, while the 
molecules are essentially equidistant in the neutral domain, so that the space and time de- 
pendent ionicity matches the staggered displacements averaged over the corresponding region. 
When the phases of the staggered displacements are almost uniform in the ionic background, 
few I-I solitons are present. Thus the neutral domain grows smoothly and rapidly (Fig. 5 
of ref. 18). Otherwise many I-I solitons are present, colliding with the neutral-ionic domain 
walls. Then, the growth of the neutral domain is rather irregular (Fig. 7 of ref. 18). In these 
examples, the coherence is strong in the former, and weak in the latter. 

In this paper, the coherence is shown to be controlled by the field strength and the pulse 
duration. When the pulse is strong and short, the ionic-to-ncutral transition is achieved with- 
out producing many I-I solitons. The charge transfer from the acceptor to donor molecules 
takes place on the same time scale with the disappearance of dimcrization. When the pulse 
is weak and long, however, many I-I solitons are produced. The dimerization-induced polar- 
ization is disordered and the staggered lattice displacements disappear on average to restore 
the inversion symmetry before the charge transfer takes place to bring the system neutral. 
Thus, the intermediate phase is regarded as paraelectric and ionic. The lifetime of the para- 
electric ionic phase becomes long as the pulse becomes weakened. In paper III, the coherence 
is shown to be indeed strong when the pulse is strong and short. The phases of the staggered 
lattice displacements remain to be almost uniform long after the photoirradiation, so that an 
interference effect is caused by a double pulse. 

In addition to the excitonic excitations summarized above, we apply pulses with infrared 
frequencies to demonstrate that the transition is induced also when the pulse resonates with 
the optical lattice vibrations. The threshold behavior is again observed. In paper III, the 
interference is shown to be closely related with the optical lattice vibrations. 

Let us come back to the coherence observed experimentally.^^) The used pulse would 
correspond to the strong and short pulse in the theoretical context. Then, the mainly ionic 
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state before the transition would have very low density of I-I solitons, some of which are 
already present by thermal fluctuations before the photoirradiation. Thus three-dimensional 
ferroelectric order initially survives for a while. After a significant density of neutral domains 
are created by the photoirradiation, if the density of I-I solitons would remain quite low 
during the proliferation of neutral domains, the phases of the neutral-ionic domain walls 
can in principle be matched (i.e., their motion is coherent), allowing to be experimentally 
observed as oscillating photoreflectance. If high density of I-I solitons were produced after the 
photoirradiation, on the other hand, the motion of the neutral-ionic domain walls would be 
irregular, making the experimental observation hard. These arguments are of course based 
on the numerically derived assumption that the motion of the neutral-ionic domain walls is 
sensitively affected by collisions with I-I solitons. 

We compare the motion of the neutral- ionic domain walls with the so-called domino effect. 
A spontaneous growth after some trigger is often described as the domino effect, but it needs 
caution. Propagation of structural changes is theoretically realized in a one-dimensional lattice 
model coupled with localized electrons after an electron is excited if the strength and range 
of the elastic coupling are appropriate.^^-* It proceeds site by site, i.e., a displacement at a 
certain site causes another displacement at the next site, then at its next site, and so on. At 
each site, the lattice is displaced under friction on the adiabatic potential. In the ionic-to- 
neutral transition, the one-dimensional dynamics is similar to this domino-like behavior only 
in the sense that, once triggered properly, a new domain spontaneously grows accompanying 
the structural changes. However, the friction is so weak that the phases of the staggered 
lattice displacements are indeed important in the transition dynamics, as shown in paper III 
in the context of the interference effect. Thus, the coherent dynamics of the present model is 
very different from the domino- like propagation of structural changes in the model of ref. 33. 
Another difference is in the itineracy of electrons in the present model, where both electrons 
and lattice displacements show cooperative and coupled dynamics. 

When comparing the simulated dynamics with the experimental findings in more detail, 
we would need to take account of three-dimensionality and energy dissipation. If inter-chain 
couplings are significant, small neutral domains in the ionic background are energetically 
unfavorable. Then sufficiently large density of neutral domains would be needed to suppress 
the effect of inter-chain couplings and to further increase the density of neutral domains. 
In the one-dimensional model used in this paper, once metastable neutral domains appear, 
they easily grow without energy loss from surroundings if the neutral and ionic phases axe 
almost degenerate in equilibrium. Therefore, we cannot quantitatively compare the threshold 
behavior. Nevertheless, neutral domains grow first along the chains, so that the present study 
would be an important step toward more detailed descriptions. 

If energy dissipation is not negligible, a part of the energy supplied by the photoexcitation 
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flows into electrons in different orbitals or different lattiec vibrations (including molecular 
ones) from those which are considered in the present model. As the transition takes more 
time (i.e., as the pulse becomes weaker), the effect of the dissipation would be larger and 
negative as long as the ionic-to-neutral transition is concerned. Indeed, the relation between 
the strength and the duration of the pulse at the threshold in the present isolated system 
deviates for weak and long pulses from that derived for a classical model with the help of a 
stochastic equation.^) It implies that the supplied energy is effectively used by prohibiting the 
energy flow into surroundings. 

In general, these effects of three-dimensionality and energy dissipation would become 
substantial for the long-time behavior. However, the recently observed short-time behavior 
would not suffer from them so much. The present result concerning the strongly nonlinear 
relation between the final ionicity and the density of absorbed photons would be qualitatively 
unchanged, although three-dimensionality or energy dissipation may quantitatively modify 
the relation. It would be worth experimental investigation to change the field strength to seek 
a possibility for changing the ionicity and staggered lattice displacements simultaneously or 
separately in a controlled manner. 
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